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We solve the Regge- Wheeler equation for axial perturbations of the Schwarzschild 
metric in the black hole interior in terms of Heun's functions and give a description 
of the spectrum and the eigenfunctions of the interior problem. The phenomenon 
of attraction and repulsion of the discrete eigenvalues of gravitational waves is dis- 
covered. 
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1 Introduction 



The perturbations of the gravitational field of black holes (BH) is a very active and 
large area of analytical, numerical, experimental and astrophysical research. Ongoing 
and future experiments, based on perturbation analysis and numerical calculations, are 
expected to give critical tests of the existing theories of gravity pi EI- For the experimental 
study of the gravitational waves a unprecedented precision of the already existing and 
future detectors must be reached PlEj- Since this is an extremely arduous and expensive 
task, an important goal of the theory is to reach detailed description of all features of the 
studied phenomena. 

For the study of the perturbations outside the BH horizon one uses the quasi-normal 
modes (QNM) with well known complex spectra QNM do not form a complete basis of 
functions in the space of the linear metric's perturbations. Nevertheless, they are widely 
used as a tool for theoretical study of the perturbations in the outer domain of BH. The 
QNM are believed to be relevant to the quantum gravity, too ^4^. 

Analogous tool for the study of perturbations of the black hole interior is not known 
even in the simplest case of Schwarzschild black holes (SBH). Since the BH-spacetimes are 
considered to describe uniform physical objects with their exterior and interior, separated 
by an event horizon, it is natural to study the perturbations of the gravitational field in 
both domains, at least by reason of completeness. 
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In the classical theory of BH it seems reasonable to exclude from the most of the 
considerations the inner domain, because it does not influence directly the observable 
outer domain. Nevertheless, the excitations of degrees of freedom in the inner domain 
may be essential for the consideration of the total BH entropy and for construction of 
the quantum theory of BH. The interior excitations may yield observable effects in the 
exterior of the BH, due to the quantum effects, see, for example, the recent article ^ and 
references therein. To discuss such problems one needs a proper mathematical description 
of the BH-interior's perturbations. 

Another reason to study the BH interior was discovered in the very recent article 
[Hj: It was shown that without excision of the inner domain around the singularity one 
can improve dramatically the long-term stability of the numerical calculations. Although 
the region of the space-time that is causally disconnected can be ignored for signals and 
perturbations traveling at physical speeds, numerical signals, such as gauge waves or 
constraint violations, may travel at velocities larger then that of light and thus leave the 
physically disconnected region. 

Thus the study of the perturbations of the BH interior becomes an important issue 
both from analytical and numerical point of view. 

In the present article we are starting the analytical and numerical investigation of 
the mathematical problem, introducing a large set of solutions for description of the 
perturbation of the SBH interior and studying the basic properties of the corresponding 
functions. We believe that this will help a more deep understanding of the problems of 
BH physics. 

We show that one can find some unexpected new results for the perturbations of the 
SBH interior. In particular, we have discovered that variations of the boundary conditions 
at the event horizon r = 1 produce attraction and repulsion of the discrete spectral levels 
of the gravitational waves in the SBH interior. 



2 General Consideration of the Interior Problem 

The proper mathematical ground for the investigation of the problem is the Regge- Wheeler 
equation (RWE) 



Its study has a long history as well as important and significant achievements [HI E] . 

The RWE describes the axial perturbations of the Schwarzschild metric of mass M in 
the first order linear approximation. In Eq. ()2.1|) x = r + r,,^^ In (|r/rg^^ — l|) denotes the 
Regge- Wheeler "tortoise" coordinate, r^^^ = 2M is the Schwarzschild radius and ^s,iit, r) 
is the radial function for perturbations of spin s and angular momentum / > s. Hereafter 
we use units rg^^=2M = l. The effective potential in RWE reads 



For r G [0, 1] the area radius r in it can be expressed explicitly as a function of the variable 
X using Lambert W-function |H]: r = W{—e^~^) + 1. 

The most important case s = 2 describes the gravitational waves and is our main 
subject here. 



(2.1) 




The standard ansatz ^sj{t, r) = i?£_s^/(r)e*^* brings us to the stationary problem in the 
outer domain r > 1: 



dlR,,,,i + {e^ - Re,s,i = 0. (2.2) 

Its exact solutions were described recently ^ in terms of the confluent Heun's functions 
HeunC{a, P,'y,6,ri,r) [TU] . 

Here we introduce for the first time the exact solutions of the Regge- Wheeler equation 
in the interior of the SBH. It is well known jT^ that because of the change of the signs 
of the components gu = —l/grr = 1 — l/r of the metric, in this domain the former 
Schwarzschild-time variable t plays the role of radial space variable r^n = t E (0, oo) 
and the area radius r G (1,0), plays the role of time variable. As a result, the Regge- 
Wheeler "tortoise" coordinate x G (— oo,0) presents a specific time coordinate in the 
inner domain. For the study of the solutions of Eq. 1)2.11) in this domain it is useful to 
stretch the last interval to the standard one by a further change of the time- variable: 

X^tin = X— l/x & { — OO, oo). 

These notes are important for the physical interpretation of the mathematical results. 
In particular, it is natural to write down the interior solutions of the RWE in the form: 

^%^i{rin,tin) = e''''^"R,,s,i {r{tin)), where 



r{Un) = W\-exp\'-^-l-jr-^\ +1 +1. (2.3) 




The dependence of this solution on the interior radial variable rj„ is simple. Its depen- 
dence on the interior time tin is governed by the RWE ()2.2|) with interior-time-dependent 
potential Vg^i- In despite of this unusual feature of the solutions in the SBH interior, this 
way we obtain a basis of functions, which are suitable for the study of the corresponding 
linear perturbations. 

Without any additional conditions the differential equation ()2.2|) has a reach variety of 
solutions which may describe different physical problems. For a correct formulation of a 
given physical problem one must restrict the class of the solutions using proper boundary 
conditions. Since the physics at the boundary influences strongly the solutions of the 
problem, the boundary conditions are a necessary ingredient of the theory and define the 
physical problem under consideration. 

The solutions of the Eq. ()2.2j) . which are not regular at the point r = are not 
stable, since they grow up infinitely as tin +oo. In addition they have a complicated 
series expansion with logarithmic terms JU]- Hence, these solutions are not single- valued 
functions and have no suitable properties for physical applications^. 

Solutions of the Eq. ()2.2p . which grow up infinitely with respect to the future direction 
of the exterior time t exist in the exterior domain of SBH, too. They are physically 
unacceptable and therefore are excluded by formulation of a proper boundary problem 
E] . Only after that the physical problem is fixed and one can prove the stability of 
the remaining solutions (QNM) and small perturbations of general form with respect to 

^It is well known that analogous solutions exist in the much more simple case of confluent hyperge- 
ometric equation. For example, because of the same reasons these are not allowed as solutions of the 
hydrogen atom problem in quantum mechanics. 



the exterior-time evolution All these solutions unavoidably grow up infinitely both 
at the space infinity and at the horizon El El • 

By analogy with the consideration in the outer domain of SBH, in the inner one we 
are considering only the solutions, which are stable with respect to the future direction of 
the interior time. Indeed, the common physical sense states that only the stable solutions 
are of physical interest. Therefore in this article we study the solutions, which are regular 
at the point r = 0, formulating the corresponding boundary problem. The remaining 
solutions are stable in the future direction of the interior time: tj„ — > +00 and have a 
nontrivial spectrum with novel properties. 

From a pure mathematical point of view in the present article we are considering the 
two-singular-points-boundary problem ^13 [Hj for the RWE on the interval r G [0, 1]. In 
this sense here we present an extension of the article , where only exact solutions of 
different types on the intervals r G [0, 00) and r G [r*, 00) C [0, 00) were studied in detail. 

There exist one more possibility: to study two-singular-points-boundary problem for 
the Eq. ()2.2|1 on the interval r G [0, cx)) using some contour in the complex plane 
with singular ends and 00, which turns round on the third singular point r = 1. Such 
solutions may describe mutual perturbations both of SBH exterior and interior. This 
interesting possibility seems to correspond to the perturbations of the Kruscal-Szekeres 
extension of the Schwarzschild solution JT]. It is beyond the scope of the present article 
and will be considered somewhere else. 

3 The Exact Solutions in the SBH interior 

The RWE ()2.2|) has three singular points: the origin r = 0, the horizon r = 1 and the 
infinite point r = oo. The first two are regular singular points and can be treated on equal 
terms. The last one is an irregular singular point, obtained by confiuence of two regular 
singular points in the general Heun's equation. As a result, in the interior of the SBH 
around the regular singular points we have the following local analytical solutions [9j^: 

i?(°)-^(r) = r^+ie^"(i°(^-")+")ifeunC(2ze, 2s, 2ie, 2e^, - /(/ + 1), r) , (3.1) 

i?«±(r) =r'^^(^''^^<^^^''^HeunC{±2iE, ±2ie, 2s, -26^, 2e^+s^-l{l + l), 1-r). (3.2) 

The solution ()3.1|) is the only one, which is regular at the point r = 0. In addition it 
is the only single- valued local solution, defined around this point. 

The solution R^l^^iir) in (j3.2|) can be derived from the one, used in jOj, applying proper 

transformations of the Heun's functions |10|. For e ^ the functions R^l^^i{r) (j3.2j) define 
a pair of linearly independent solutions with the symmetry property 

R%ir)=R%^,{r). (3.3) 
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Figure 1: Two solutions of the continuous spectrum. 

4 The Interior Spectral Problem 
4.1 The Continuous Spectrum (e^ > 0). 

For real > the problem has a real continuous spectrum. The corresponding solutions 
Re,2,2{tin) are illustrated in Fig. [T] They define a basis of stable normal modes for 
perturbations of the metric in the SBH interior and can be used for Fourier expansion of 
perturbations of more general form. 

As seen in the Fig. the use of the inner time variable tin instead of the tortoise 
coordinate x makes the variations of the solutions almost equidistant in the limit tj„ 
—oo. This property makes the inner time variable tin to seem the most natural one in 
the SBH interior. Using variables like r or x we lose this property and meet numerical 
problems in this limit. 



4.2 The Special Solutions with = 0. 

The zero eigenvalue e = yields a degenerate case in which the two functions ()3.2|) 
coincide. Then we obtain the infinite series of polynomial solutions: 

R%s,i{r) = p^,;,(|~2i' ^_^^ r^'JacoUP{l-s, 2s, 0, l-2r) = r' + . . . , (4.1) 

which are finite on the whole interval r G [0, 1]. These functions form an orthogonal basis 
with respect to the measure w{r) = (1 — r)^*. Here JacobiP{n,a,b, z) stands for the 
standard Jacobi polynomials and Poch{c, d) - for the Pochhammer symbol. 

The basis of the stable normal modes 1)4.11) is suitable for series expansion of the 
perturbations of e = 0. These solutions do not depend on the variable r^^ and are the 
only ones, which are finite at both singular ends of the interval r G [0, 1]. 



4.3 The Discrete Spectrum (e^ < 0). 

The discrete spectrum of the problem at hand corresponds to pure imaginary values of 
the parameter e, i.e. to Laplace transform of the perturbations of general form, instead 
of their Fourier transform, which corresponds to real values of e. Thus, studying the the 

^For a detailed explanation of the notations and physical meaning of the solutions one can consult 
this article 



case < we are constructing a basis for Laplace expansion of the perturbations of more 
general form in the SBH interior. 

For a qualitative analysis of the discrete spectrum we can apply the Ferrari-Mashhoon 
transformation JH]: x ^ x = —ix, fg = —ivg, V ^ V = —V. It reduces our 

spectral problem to a standard one: to find the bound states in the inverted potential 
V. For s > 2 the latter describes two negative-valued potential wells i) a finite one, in 
which one may have several negative eigenvalues, and ii) an infinite one, in which one has 
infinite series of negative eigenvalues. In both cases < 0. Hence, Re (e) =0. Then the 
functions ()3.1|) . ()3.2|) and all their parameters are real. In particular. 



(1)-, 



(4.2) 



~ cos a R^l'^^iir) + sin a R^'J^^^iir) 
with real transition coefficients r'j'!^ 9] and real mixing angle a: 

tana = r?l(£:, s,/)/r'j'^(£:, s,/) G M, for e = iei. 



(1)- 



(4.3) 



The physical meaning of the mixing angle a is obvious: since the solutions R^^^^i{r) 



describe the ingoing into the horizon waves and the solutions R 



(1)- 

£,S,l ^ 



- the outgoing from 



the horizon waves [9j, this angle describes the ratio of the amplitudes of these waves in 
their mixture ()4.2j) . The value a = corresponds to absence of the outgoing waves, and 
the value a = tt/2 describes the case without ingoing waves. Thus, fixing the value of the 
mixing angle a we are defining completely the physical problem under consideration. 

When r — >■ 0, the asymptotic of the inverted potential is ~ (1 — s^)r~^. Hence, we 
arrive at the well known from quantum mechanics specific physical problem of "falling 
at the centrum" - the singular point r = 0. Mathematically this means that the corre- 
sponding differential operator on the interval r G [0, 1] has a defect ^3] and we need an 
additional physical conditions to fix the eigenvalue problem. In our case such conditions 
are the regularity condition of solutions at the singular point r = and the fixation of 
the mixing angle a, i.e. the justification of the physical conditions at the second singular 
point - the horizon r = 1. Thus, for any given value of a the Eq. ()4.3j) presents a specific 
spectral equation for the eigenvalues. 

Unfortunately, at present the transition coefficients r'j'^(£:, s, I) are not known explicitly 
[TUj . Therefore we replace the spectral equation ()4.3|) with the equivalent one: 



(^1 ) cos a (ri ) + sin a R\, [n , 
R%{t^) cos a ^%{r^)^^v^ a R^l^ir,) 



0, 



(4.4) 



where ri 2 G (0, l),ri 7^ r2 are two otherwise arbitrary points. Because of the symmetry 
property ()3.3|) . it is enough to study only the solutions en,s,i of the equation ()4.4|) with 
positive Ej = Im (£«,«,/) > 0. 

We apply the complex plot abilities of Maple 10 to find a relatively small vicinity of 
the complex roots of this equation and then we justify their values by Muller's method 
for finding of complex roots Our numerical results for Im (en,2,2) are shown in the 

Table □ and in Figs. EHi 



In the Fig. 2 we see the first 18 eigenvalues (including the zero eigenvalue) of the RWE 
for SBH interior for s = I = 2 and a = 0. The two series: n = 0, 1, ... 6; and n = 7, . . . are 
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Im (£:„,2,2) 


n 




n 




1 


0.7523 


7 


4.5138 


13 


7.5657 


2 


1.3463 


8 


5.0256 


14 


8.0714 


3 


1.8999 


9 


5.5357 


15 


8.5768 


4 


2.4361 


10 


6.0446 


16 


9.0823 


5 


2.9627 


11 


6.5524 


17 


9.5867 


6 


3.4833 


12 


7.0594 


18 


10.0939 



Table 1: The values of the first 18 eigenvalues of the RWE for SBH interior for s = I = 2 
and a = 0. 




Figure 2: The first 18 eigenvalues (for s = I = 2 and a = 0). 

clearly seen. They are placed around the straight lines ej = n/2 + 1/2 and ej = n/2 + 1, 
corresp ondingly. 

In FigElwe see a specific attraction and repulsion of the eigenvalues en,s,i{c()- Such 
behavior is typical for the eigenvalue problems related with the Heun's equation JUI- It 
is well known in the quantum physics. To the best of our knowledge, this phenomenon is 
observed for the first time in the theory of SBH. Several additional remarks are in order: 

1. Solving Eq. ()4.4|1 one is unable to obtain the eigenvalues £Ar/2,2,2(f ) = ^^^^ ^ '^'^ ^ 
i.e., when a = |/ mod vr (the small circles in Fig|^. A pure numerical study of the 
RWE ()2.2j) shows that in these cases the solutions have the same behavior as the one, 
demonstated in FigOl Hence, one has to be careful, taking the corresponding limits in 
the analytical expressions (jSUl), ()3.2p . 

2. In the second series: n = 7,... the dependence of the eigenvalues on the mixing 
angle is similar to the one, shown in FigEJ but Im (£,1,2,2) are decreasing functions of a. 
The amplitudes of their variations are extremely small (< 10^^) and they are practically 
constant on the larger part of the interval a G [0, 2tt]. 
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Figure 3: From right to the left: The first 18 eigenfunctions (for n=l, 2,. . . ,18, s=l=2 
and a = 0). 
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Figure 4: The dependence of Im (£0,2,2) on the mixing angle a e [0, 27r]. 



3. There exist two more series of eigenvalues en,2,2- For them all Im (£n,2,2) are 
negative. These have the same absolute values as the eigenvalues, shown in Table ^ and 
in FigEJ Their dependence on the mixing angle a is shifted by ^ and correspondingly 
inverted. Because of this shift, the behavior of the corresponding eigenfunctions is the 
same as the one, demonstrated in FiglHl This reflects the symmetry property ()3.3j) of the 
solutions. 

4. The behavior of the eigenfunctions -Ri,2,/('") with 1 = 3,... is the same as of the 
functions -Ri,2,i('^) in FigOl The dependence of the eigenvalues £1,2,; ('^ = 0) on the angular 
momentum / is shown in FiglSJ 

5. The eigenfunctions Rn,2,i{f) with n > 1,1 > 2 start from infinity at r = 1 and 
approach zero at r = 0, oscillating around the zero value. The number of oscillations 
depends on the numbers n and /. 
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Figure 5: The dependence of the eigenvalues ei,2,«('^ — 0) on the angular momentum /. 

6. Our results bring to light a complementary information about the stability problem 
of SBH. Up to now only the stability of the exterior domain has been studied in linear 
approximation ^2)- Here we have found the following details about stability properties 
of the solutions in the SBH interior: 

• As we see in Figs. ^ |Sl all solutions of Eq. (j2.H) . which are considered in the 
present article, are stable in direction of the interior-time-future (t^ ^ +00), according 
to Lyapunov criterion. If e ^ 0, these solutions are unstable in direction of the interior- 
time-past (tin —00) according to Lyapunov criterion. 

• There is a basic difference in the stability properties of the solutions of continuous 
spectrum (real e, > 0), which are of limited variation when tin ~oo and thus - 
stable in direction of the interior-time past in this weaker sense, and the solutions of 
discrete spectrum (imaginary e, e"^ < 0), which grow up infinitely when tj„ — —00 and, 
hence, are unstable in this direction of the interior time in any sense. 

• All other interior solutions of the Regge- Wheeler equation with e ^ 0, which are 
not considered in this article, are unstable in the future direction of the interior time and 
drop out of the physical consideration. 

• There do not exist interior solutions, different from the above ones with e = (see 
Eq. fl4.1|) ). which are stable in both directions of the interior time. 



5 Concluding Remarks 



In this article we have shown that the Heun's functions are an adequate powerful tool for 
description of the linear perturbations of the gravitational field not only outside the SBH 
horizon but in their interior domain as well. 

It is well known that every analytical function is defined completely by its singularities 
in the whole compactified complex domain. In the previous investigations of the specific 
properties of the solutions of RWE was used only the structure of these singular points, 
see the review articles [3J and the huge amount of reverences therein. Using many different 
mathematical technics, this way were reached important results for perturbations in the 
exterior of SBH. Some of the used methods of that type may be not applicable for the 
study of perturbations in the SBH interior. 

Here we overcome the limits of that methods, using directly the properties of Heun's 
functions, see for example Eq. ()3.1|) . ()3.2|) . ()3.3|) . ()4.3|1 . etc. and writing down the corre- 
sponding explicit analytical solutions of the RWE. 

All numerical results were obtained using Maple 10 package for performing direct 
calculations with these functions. Owing to the new version 10 of the Maple package, 
at present we are able to work with the Heun's functions in the complex domain, thus 
obtaining a new freedom and new tools for the study of different new boundary problems 
for RWE, as well as for study of the old ones, applying new technics. 

Using Heun's functions we have obtained the analytical and the numerical description 
of all solutions of the Regge- Wheeler equation in the SBH interior, which are stable 
in direction of the interior-time future and single-valued functions in the whole complex 
domain. Having in our disposal these solutions we are ready for a more deep understanding 
of the problems of BH physics. 

For example, as we have seen, the perturbations of SBH metric do not change the 
singularity at the point r = 0, at least in the first order of perturbation theory. This result 
can be used for further justification of the numerical treatment of BH problems without 
excision of the interior domain. Now it becomes clear, that due to the unavoidable 
numerical errors in the interior of the BH one will lose the small stable single-valued 
physical solution around the point r = at the background of the big nonphysical solution. 
Hence, in order to obtain physically correct and stable numerical results one must suppress 
by proper techniques this nonphysical solution. 

We have shown that under proper boundary conditions there exist interior solutions 
both of real continuous spectrum and of pure imaginary discrete spectrum. The first 
may form a basis for Fourier transform of linear perturbations of general form in SBH 
interior. The second may form analogous basis for Laplace transform. The investigation 
of the completeness of these sets of functions as a basis for the corresponding expansions 
and the analytical proof of the numerical results of the present article remain an open 
problem. 

Our numerical study shows that there do not exist other solutions of the spectral 
equation (j4.4j) . In particular, in the two-singular-points-boundary problem on the interval 
r G [0, 1] we do not found solutions with Re{e) x Im{e) ^ 0, in contrast to the case of 
exterior domain perturbations of SBH. 

The discrete pure imaginary eigenvalues have a nontrivial behavior under changes 
of the boundary conditions at the horizon, described by the mixing angle a. The phe- 
nomenon of attraction and repulsion of these eigenvalues has been observed. 



The basic properties of the considered solutions have been studied and a new infor- 
mation about stabihty of the solutions of RWE in the SBH interior was obtained. 

In the present article we have studied only the most important from physical point of 
view perturbations of SBH with s = 2. The detailed study of the solutions with s = 0, 1 
may be useful for numerical relativity in the spirit of the article jHj. 
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